Optimization of soliton ratchets in inhomogeneous sine-Gordon systems 
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Unidirectional motion of solitons can take place, although the applied force has zero average 
in time, when the spatial symmetry is broken by introducing a potential V{x), which consists 
of periodically repeated cells with each cell containing an asymmetric array of strongly localized 
inhomogeneities at positions Xi. A collective coordinate approach shows that the positions, heights 
and widths of the inhomogeneities (in that order) are the crucial parameters so as to obtain an 
optimal effective potential Uopt that yields a maximal average soliton velocity. Uopt essentially 
exhibits two features: double peaks consisting of a positive and a negative peak, and long flat 
regions between the double peaks. Such a potential can be obtained by choosing inhomogeneities 
with opposite signs (e.g., microresistors and microshorts in the case of long Josephson junctions) that 
are positioned close to each other, while the distance between each peak pair is rather large. These 
results of the collective variables theory are confirmed by full simulations for the inhomogeneous 
sine-Gordon system. 

PACS numbers: 05.45.Yv, 85.25.Cp, 73.40.Ei 
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I. INTRODUCTION 

Ratchet or rectification phenomena appear in 
many different fields ranging from nanodevices 
to molecular biology [H S S S| • In the sim- 
plest model a point-like particle is considered 
which is driven by deterministic or non-white 
stochastic forces. Under certain conditions re- 
lated to the breaking of symmetries, unidirec- 
tional motion of the particle can take place al- 
though the applied force has zero average in 
time. 

These particle ratchets have been generalized 
to spatially extended nonlinear systems, in which 
solitons play a similar role to the above point 
particles 0, H, H, HO, U) I^ particular, 
solitons in nonlinear Klein-Gordon systems have 
been shown to move on the average in one di- 
rection, although the driving force has zero time 
average, if either a temporal or a spatial symme- 
try is broken. 

In the first biharmonic driving force has 

been used which breaks a time shift symmetry 
TsL l3 | . Here the mechanism of the ratchet effect 



has been clarified in detail by a collective vari- 



able theory 0, 16||, which uses the soliton po- 
sition X{t) and width l(t). Due to the coupling 
between the translational and internal degrees 
of freedom, energy is pumped inhomogeneously 
into the system, generating a directional mo- 
tion. The breaking of the time shift symmetry 
gives rise to a resonance mechanism that takes 
place whenever the width l{t) oscillates with at 
least one frequency of the external ac force. This 
ratchet effect has been confirmed by experiments 
with annular Josephson junctions 17[ which can 
be modeled by sine-Gordon systems; here flux 
quanta (fluxons) play the role of the solitons. 
As the external ac force, biharmonic microwaves 
have been used. 

The biharmonic force has recently been re- 
placed by more general periodic forces which 
can be expressed by Jacobian elliptic functions. 
There, it turns out that the average soliton ve- 
locity exhibits an extremum for a certain value 
of the modulus of the elliptic functions . This 
means that there is an optimal form for the pe- 
riodic force. 

Another way to obtain a soliton ratchet is 
to break the spatial symmetry. This has been 
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demonstrated very recently by the introduc- 
tion of point-like inhomogeneities into nonlinear 
Klein-Gordon systems [19. , 2D.] ■ These inhomo- 
geneities can be modeled by (5-functions, if their 
spatial extent is much smaller than the char- 
acteristic length for the system (the Josephson 
penetration length in the case of the Josephson 
junctions). In order to achieve a ratchet effect 
the (5-functions must form asymmetric periodic 
arrays. This asymmetry translates into an asym- 
metry of an effective potential U{X, I) which ap- 
pears in a collective variable theory with the 
variables X{t) and l{t). 

As an example, three (5-functions of equal 
strength have been positioned in each cell of an 
array with period L. The above theory yielded 
a nearly perfect agreement with the simulations 
for the inhomogeneous sine- Gordon system [2fl| . 
By choosing different arrangements for the po- 
sitions of the (5-functions within the cells of the 
array it was possible to increase the average soli- 
ton velocity (v) which means that the ratchet 
system provides better transport. 

The aim of this paper is to search for the 
optimal shape and arrangement of localized in- 
homogeneities, in the sense that (v) becomes 
as large as possible. This means that in- 
stead of (5-functions other functions which repre- 
sent localized inhomogeneities (e.g., Gaussians, 
Lorentzians, box-type functions etc.) have to be 
tested, and that the arrangement of these inho- 
mogeneities within the cells, but also the length 
L of the cells in the periodic array, all have to 
be optimized. 

In section II nonlinear Klein-Gordon systems 
with localized inhomogeneities are introduced. 
In section III a collective coordinate (CC) theory 
is developed which results in a set of ODEs con- 
taining an effective potential U which depends 
on the inhomogeneities. It turns out that their 
most significant features for maximizing (v) are 
their positions, heights and widths, while their 
detailed shape is unimportant. For this reason 
we work with the simplest shape, namely that of 
boxes, where U can be calculated analytically. 
In section IV the optimal effective potential Uopt 
is obtained by using an ansatz in terms of Ja- 
cobi elliptic functions. Uopt essentially exhibits 



two features: double peaks consisting of a pos- 
itive and a negative peak with a steep slope in 
between, and long fiat regions between the dou- 
ble peaks. In section V we show that the essen- 
tial features of Uopt can be produced by choosing 
appropriately the positions, heights and widths 
of box inhomogeneities. In this way very high 
soliton velocities can be obtained. These results 
of the CC- theory are confirmed by full simula- 
tions for the sine-Gordon system with the above 
inhomogeneities (section VI). 

In the last section we stress that an impor- 
tant step in the optimization consists in using 
inhomogeneities with opposite signs. Both types 
of inhomogeneities were already used in long 
Josephson junctions, namely microresistors 
(critical current X decreased) and microshorts 
(Jc increased) 2X1- However, so far only 
equidistant arrays of either type were used 
in experiments. We propose to use a specific 
combination of both types, resulting in an asym- 
metric periodic array, which yields very high 
average soliton velocities both in a collective 
coordinate theory and in full simulations for the 
inhomogeneous sine-Gordon system. 



II. INHOMOGENEOUS NONLINEAR 
KLEIN-GORDON SYSTEMS 



A ratchet effect has recently been obtained 
2ol | by breaking the spatial symmetry of 
nonlinear Klein-Gordon systems with inhomo- 
geneities introduced via a potential V{x) 



+ 



+ 



dU 



[1 + V{x)]=f{t). 



Here 0(x, t) is a scalar field, (px and (pt are par- 
tial derivatives with respect to space and time, /? 
is a damping coefficient, U{4>) = 1_— cos0 is the 
sine-Gordon (sG) potential , and U = |(0^ — 1)^ 
the 0^-potential. As both cases have produced 
very similar results, Ecl] we will concentrate on 
the sG-model in the following, but we will always 
indicate which results also hold for the (/)^-model. 
f{t) = Asm{ujt + 6o) is an external ac force with 
amplitude A, frequency lu and initial phase Sq. 
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V{x) consists of periodically repeated cells of 
length L; each cell n contains an asymmetric ar- 
ray of strongly localized inhomogeneities gi{x) 
which are placed at positions Xi within the cell; 
i.e. 



Vix) = 9i{x -Xi- nL). 



(2) 



So far only the case of point-like inhomo- 
geneities was considered which can be modeled 
by 5-functions. Taking three 5-functions with 
equal strengths and choosing their positions Xi 
in certain asymmetric ways, a sine-Gordon soli- 
ton (also named kink in the following) moves 
either to the right or to the left on the aver- 
age 01 . Depending on the choice of the driving 
frequency qualitatively different results were 
obtained for the modulus v of the average soli- 
ton speed (w) as a function of the driving am- 
plitude A: a) For the low frequency uj = 0.015 
there is an up and down staircase with the maxi- 
mum average velocity Vmax = 0.03820, b) for the 
intermediate frequency u = 0.05 there are five 
"windows" (regions of A with constant v > ) 
with fmax=0. 03183, c) for the relatively high fre- 
quency UJ = 0.1 there is only one window with 
height Vmax = 0.063666. For even higher fre- 
quencies, V vanishes for all amplitudes A. The 
results for (p^ nonlinearity are quite similar. 

Vmax is always small compared to the critical 
velocity c which is unity for the dimensionless 
system (|T)). For long Josephson junctions c is 
the Swihart velocity 2^ . 

The aim of this paper is to optimize both the 
shape and the array of the inhomogeneities in 
the sense that v becomes as large as possible. 
Naturally, there are still other criteria for the 
optimization. E.g., one can try to maximize the 
area under the curve v vs. amplitude A ; this will 
be discussed later. Other optimization strategies 
include studying the energy efficiency. We do 
not consider this approach here, but we refer the 
reader to Ref.[23i]. 

In order to achieve any kind of optimization, it 
is necessary to have a theory which allows calcu- 
lation of (v) for a given array of inhomogeneities. 
A CC theory with two variables, namely position 



X{t) and width l[t) of the soliton, has turned 
out to be very successful in the sense that 
in the case of 5-functions the results for (v) as a 
function of the driving amplitude A agree very 
well with the simulations for the original system 
(full numerical solution of Eq. (Q)). Deviations 
between theory and simulation are found only in 
the case of strong driving {A = 0{1)). For even 
stronger driving the theory is no longer appli- 
cable, because in the simulations kink- ant ikink 
pairs appear spontaneously and scatter with the 
soliton under consideration. The 2-CC theory is 
characterized by two ODEs for X{t) and l{t), 
which contain forces that are defined as par- 
tial derivatives of an effective potential U{X,l), 
which is calculated and discussed in the next sec- 
tion. 

We take for our computations a rather large 
damping, choosing (3 = 1, for the following rea- 
sons: 

a) Our CC-theory does not take into account 
the phonons which are radiated by the soli- 
ton due to its acceleration by the driving force. 
These phonons are quickly damped out when the 
damping is large. 

b) We want to avoid several phenomena which 
can occur for very small damping: e.g., chaotic 
behavior and current reversals. 

Finally, we remark that the spatial symmetry 
can be broken also by the introduction of an ad- 
ditive inhomogeneity 7(x) in Eq^ instead of our 

riTJ 

multiplicative inhomogeneity ^^V"(a;). In the 

case of long Josephson junctions an additive in- 
homogeneity has been experimentally realized by 
a current injection with a profile 7(0:) [1^. This 
produces an effective potential which is similar 
to an asymmetric potential. 

III. COLLECTIVE COORDINATE THEORY 

For the unperturbed sine-Gordon equation, i.e. 
Eq. without damping, inhomogeneity and 
driving, the 1-soliton solution reads 



(l)o{x,t) = 4 arctan[exp(7 



X 



vt^ 



(3) 
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where v is the sohton velocity, /q = 1 is the rest 
width, and ^o/7 = ^oVl ~ is the Lorentz con- 
tracted width. 

For perturbed nonhnear Klein-Gordon systems 
the so-called Rice ansatz jl^ turned out to be 
very successful. E.g., for the ac driven 0^- 
system an unexpected resonance was predicted 
(and confirmed by simulations), which is situ- 
ated in the gap below the phonon spectrum at 
half the freq uency of the internal mode of the 0^- 
kink (26l.l27l|. The Rice ansatz for the sG-systems 
reads 



X — Xit) 

X, /) = 4 arctan{exp [—^^] } (4) 

with soliton position X{t) and width l{t). The 
2-CC theory for the perturbed sG-equation (P) 
yielded a set of two coupled ODEs 



X X XI 

Mo/o- + /3Mo/o--Mo/o-^ = F"^ + F*"' (5) 



i i x'^ 

aMok- + paMok- + Mok 



(6) 



I ' ' ' 

where a = 7r^/12, Mq = 8, /q = 1 and 

with the topological charge q = 27i. 



K 



int 



dE 



is a force which arises from the soliton energy 

E{X,l,i) = ljMoX'+l^jaMoP+lMo{^j+^J 

(9) 

including internal energy due to width oscilla- 
tions. 

There are two forces which appear due to the 
potential V{x) in Eq. ((T)), viz. 



°° dU . 



dU 
dX 



(10) 



K 



inh 



dU 

dx—-V(x) : = 
ol 



In this way the effective potential 



U{X,l) 



dxtj{4>)V{x) 



(12) 



is introduced that characterizes the influence of 
the inhomogeneities on the soliton dynamics. 
We remark that all the above results also hold for 
the 0^-model: In this case = tanh[(x — X)/l] 
and the parameters are g = 2, Mq = 2^/2/3, Iq = 
y/2, and a = (tt^ - 6)/12. 

Before we specify V{x)^we evaluate U (0) . First 
we derive the relation f/(0o) = {^y/{2-f^) = 
2sech^(7z), with z = x — vt, by inserting 0o 
into the unperturbed sG-equation. However, 
[/(0) can not be obtained in this way because 
(j){x, X, I) is not an exact solution, but rather an 
ansatz. However, we_can generalize the above re- 
lation for f/(0o) to L/(0) = i/V^ = 2sech\z/l), 
with z = X — X , and this can be verified by e.g. 
MATHEMATICA |28|. 

Thus we finally obtain 



U{X,l) 



da;2sech2^-^- V^(x). (13) 



Now we can insert the superposition (jSJ of lo- 
calized inhomogeneities Qi and obtain 



u{x,i) = Y,Y.u, 



in) 



(14) 



(n) 



axisech — - — 
I' 



gi{x -Xi- nL). 

(15) 

We can check that the effective potential is in- 
deed periodic, U{X + L,l) = U{X,l), because 
the sum goes over an infinite number of unit cells 
n of length L. 
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Our next step is to evaluate U^^'^ for vari- 
ous strongly localized functions gi and to dis- 
cuss which are the most important features of 
the inhomogeneities. We have tested Gaussians, 
Lorentzians and box functions and it turns out 
that there are only two important features of Qi, 
namely the height and width. For the above 
three cases, the U^""^ are always bell-shaped, 
and the desired asymmetry of U depends most 
strongly on the relative positions of the U-'^^ 
within the cells. In other words, the most impor- 
tant features of the inhomogeneities Qi are their 
positions Xi, then come their heights and widths, 
while their detailed shape is unimportant. 

For this reason, we can choose the simplest 
case, namely the box functions, which has the 
additional advantage that U (X, I) can be calcu- 
lated analytically (in contrast to the other two 
cases). The box function at Xi in cell n = is 
defined as 



hi for Xi-bi<x<Xi + bi 



otherwise 



This yields 



u\°'' = 2hil sinh sech 

X sech J (17) 

and we can check that in the limit hi ^ oo, bi ^ 
with finite := 2bihi the result for a 6- 
function inhomogeneity [19| can be regained: 



(0) 



2ejSech [{X — Xi)/l]. The final result 



for the effective potential (fT^ is 



^ ^ 21). 2i Z 

U{X, I) = 2 hi I sinh — ^ sech— p sech— 



with Z± = X — Xi — nL ± 6j. 

In this way the optimization problem is solved, 
but only in principle: We choose a set of box 
inhomogeneities at the positions Xi within the 
cells, numerically solve the ODEs (0), (0) for 



X{t),l{t), and compute the average soliton ve- 
locity, defined by 



{v) = {X) 



lim 



(X(t)-X(O)) 



(19) 



choosing a sufficiently long integration time. 
Then the procedure has to be repeated for many 
values of the driving amplitude A (keeping the 
driving frequency uj fixed), in order to find the 
maximal soliton speed 



Max{\{v)\} 



(20) 



for that set of parameters. Finally, we would like 
to obtain Vmax for a low, a medium and a rela- 
tively high frequency, because it is known from 
the case of the 5-function inhomogeneities that 
the curves {v) vs. A differ qualitatively for the 
above choice of frequencies (see section II). 

In practice the whole procedure cannot work 
because the dimension of the parameter space 
is far too high. E.g., if one takes three inho- 
mogeneities per cell, as in the case of the 5- 
functions, one has 10 parameters, namely three 
positions, heights and widths of the boxes, plus 
the cell length L. As one has to make sweeps 
in A in order to get Vmax for every parameter 
set, it is practically impossible to explore the 
10-dimensional parameter space. Obviously, we 
have to find a more efficient way to obtain an 
optimal set of inhomogeneities. 

IV. OPTIMAL EFFECTIVE POTENTIAL 

Our approach is not to work directly with 
the potential V{x) in the inhomogeneous sG- 
equation ((T)), but to design first an optimal 
effective potential Uopt for a given driving 
force, i.e. for fixed A and uj. When Uopt has 
been obtained, it can be represented (approxi- 
mately) by the superposition (fTH|) of the pulse- 
shaped contributions U^^\ I.e., the many pa- 
rameters in Eq. (|18|) could be determined for 
instance by a least-square fit of U to Uopt-, which 
is much more efficient than a search in the high- 
dimensional parameter space. 
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For simplicity we will not try to find Uopt{X, I) 
within the 2-CC approach. It will turn out to be 
sufficient to find an optimal potential within the 
1-CC theory which starts with the ansatz 

0o(2:,X) =4arctan{exp[7(x-X(t)]} (21) 
with 

7 = 1/Vl -X2. (22) 
This yields one ODE Q 

'y^MoX + -f/3X = F"" + F'""^ , (23) 
with F^"'' = -dU/dX and 

/oo 
dxU{(t))V{x) 
-oc 

/CO 
dx2sech^[-f{x - X)]V{x). (24) 
-oo 

Here, 'J^Mq is the so-called longitudinal relativis- 
tic mass. So far only the non-relativistic limit 
(7 = 1) of Eq. (P^j) was considered in the con- 
text of the ^-function inhomogeneities [l^ liO ]. 
Here we need the relativistic version, because 
our optimization will yield velocities which are 
no longer small compared to the critical velocity 
c = 1. 

As noted above, we do not consider here the 
expression for U in terms of V{x), but we 
want to find quite generally the optimal U for a 
given driving force F"-'^. Uopt must be asymmet- 
ric and periodic. As the most general periodic 
functions are the Jacobi elliptic functions, we use 
the ansatz 

Uopt{X) = —e ■ sn{kX, m) ■ cn{kX, m), (25) 

where m is the modulus of the elliptic functions, 
k = 2K{m)/L is a generalized wave number and 
K{m) is the complete elliptic integral of the first 
kind. 

This type of ansatz was first introduced in 
[l8| as a generalized driving force to break the 
temporal symmetry in ratchet systems, namely 
fit) = Feiiipit) = e ■ sn{flt, m) ■ cn{Qt, m) (where 



VL = 2K[m)/T with period T). Using only sym- 
metry arguments, and taking into account the 
first two terms of the Fourier series of FeiHp{t), 
it was shown in 18] that the optimal value for 
m is 0.960057, independent of the details of 
the ratchet models. However, in our case we 
will show that the situation is more complicated 
when such an ansatz is taken for the potential 
U{X) in the particle ratchet ((221) • 

The negative sign in front of e > in (j25j) is 
chosen in order to achieve a positive (f ), i.e. av- 
erage motion to the right. The specific combina- 
tion of elliptic functions is chosen because 
there is no ratchet effect in two limiting cases: 
For m = the effective potential is sinusoidal, 
and for m = 1 it is zero (except for a set of points 
with Lebesgue measure zero). Thus we can ex- 
pect that an optimal value mopt naturally exists 
in the interval < m < 1, under the condition 
that there is any ratchet effect for the chosen 
set of parameters e, L in U{X) and A, uj in the 
driving force. 

Our problem now is to find the optimal mod- 
ulus m, i.e. the optimal shape of the potential 
U{X) in the ODE (UHl), which formally is a rel- 
ativistic particle ratchet model. In Fig. ^ [v) 
vs. A is plotted for various values of m and we 
can see that the maximal speed Vmax gradually 
increases with increasing m. However, when m 
is close to one, there is a dramatic increase of 
Vmax- The reason for this can be discussed by 
considering the change of shape of U {X) in Fig. 
121 When m is increased (but is not yet close to 
1) the negative slopes in U{X) become steeper 
and the positive slopes become weaker. Both 
together increases the asymmetry of the poten- 
tial; thus the symmetry breaking is augmented, 
which explains why Vmax grows. In fact, U{X) 
looks like a (smoothed) sawtooth potential that 
has widely been used in the literature. 

However, when m becomes close to one, a new 
feature appears in U{X)^ namely long fiat re- 
gions. Their lengths increase as m gets even 
closer to one, so eventually U{X) essentially ex- 
hibits two features: 

1.) double peaks, each consisting of a positive 
and a negative peak with a steep slope in be- 
tween. 
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FIG. 1: Mean soliton velocity vs. amplitude A of the ac 
force for different moduli m of the elliptic functions in the 
ansatz Ij25|l for the effective potential. The parameters used 
are L = 4, /3 = 1, tj = 0.015 and e = -2. 
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FIG. 2: Effective potential 12511 for different values of the 
modulus m of the elliptic functions; L = 4, e = —2. 



2.) long flat regions between the double peaks. 

A potential with similar features was proposed 
for a ratchet system driven by a multiplicative 
white noise to enhance the coherent transport 
■ 

Before we proceed with our goal, namely to 
find the optimal inhomogeneities which produce 



J 



-'step ) 



i,j : integer. 



(26) 



the above features, we want to discuss the stair- 
case structure of (v) vs. A in Fig. 1. This struc- 
ture is described by 

It 

The step height Vstep = L/T does not depend on 
the shape, but only on the period of U{X), and 
on the period T = 27r/c<j of the ac force. 

Formula ()26|) should hold quite generally ac- 
cording to most of the literature ( |J| and refer- 
ences therein). However, Eq. (j26|) was explic- 
itly calculated only for a special case, namely 
an asymmetric sawtooth potential and a piece- 
wise constant driving force js^ . Our results con- 
firm Eq. (j2ni), and so do the results for the ef- 
fective potential stemming from 5-function inho- 
mogeneities jl^. Moreover, we observe that {v) 
changes in steps also as a function of other pa- 
rameters. For the following it will be important 
that (f ) as a function of m first increases, then 
decreases in steps (Fig. This means that 

Vmax is reached for a very narrow range of m, 
which we denote by [mop^]. The larger A, the 
closer [niopt] is to 1. For a larger value of the 
frequency, u = 0.05, the same effect shows up, 
namely it is found that the ratchet transport is 
enhanced for m very close to 1 provided that A 
is not too small. 

V. OPTIMAL INHOMOGENEITIES 

Our approach was to construct an array of in- 
homogeneities which produces an effective po- 
tential whose shape is approximately equivalent 
to the shape of Uopt, i-e. it must exhibit the two 
essential features found in Uopt in the previous 
section. The required effective potential can be 
obtained in four steps (here we skip working with 
U{X) from the 1-CC theory and proceed directly 
with U{X,l) in Eq.(P|)): 

1.) A single box inhomogeneity produces the 
single peak (fTTj) . two boxes with equal width 
(6i = 1)2 = h) and equal height, but opposite 
signs {hi = h,h2 = —h), produce the desired 
double peak structure, if they are placed at po- 
sitions xi, X2 not too far separated. 
2.) As to the second feature of Uopt, the long fiat 
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FIG. 3: Mean soliton velocity vs. the modulus m of the 
elliptic functions in the ansatz 12511 for the effective potential; 
L = 4, A = 0.6, 13=1, L0 = 0.015, e = -2. 



regions between the double peaks can easily be 
obtained by choosing a large period L, compared 
to the width of the double peaks. For the widths 
discussed below, a good range is L = 4 — 8; here 
(f ) depends only weakly on L. 

3. ) We need to make the negative slope within 
a double peak very steep, and the positive slope 
between double peaks very weak, because these 
are features of Uopt- Both features can be 
achieved by making the two boxes narrower and 
higher (keeping the area fixed). Unfortunately 
the effect of this procedure is limited by the fact 
that the two peaks u[^^ and 1/2^^ which are pro- 
duced by the two boxes, have a minimum width 
of about 1 (for / ~ 0(1)), see Eq. (HH). This 
width is practically the same as that of the sech"^- 
peak produced by a (5-f unction inhomogeneity, 
see below Eq. (ITTj) . 

4. ) Because of the limitations of step 3, we po- 
sition the two boxes closer together. This would 
make the negative slope steeper, if there were not 
a partial compensation of the positive and neg- 
ative contributions u'f'' and u!^^ which reduces 
the heights of the two peaks. One can compen- 
sate this effect by further increasing the height h 
of the boxes and thus the peak heights. The op- 
timal double-peak structure is produced by two 



very narrow boxes which are as close together 
as possible, namely touching each other without 
overlapping (e.g., Xi = 0.8, 0:2 = 1.0,6 = 0.1). 
Then the height can be increased as long as Vmax 
is reached for A < 1 (for the above example 
h = 25, see Fig. Ej). There is no value in in- 
creasing the amplitude of the driving force be- 
yond A ~ 1 because in the simulations for the 
original PDE kink-antikink pairs may then 
appear spontaneously and degrade the ratchet 
effect (see also end of section II). 

The above four steps together yield an enor- 
mous gain for Vmax in the order of 300%: In 
Fig. 13^ the maximal value of (v) is 0.1528 for 
uj = 0.015 , which has to be compared with the 
best value 0.0382 for the case of three (5-function 
inhomogeneities with the same uj . The value 
0.1528 for Vmax is indeed quite high taking into 
account that this is an average over positive and 
negative velocities, cf. Eq. ()19j) . Even larger val- 
ues of Vmax can be obtained by choosing a small 
damping parameter, but we do not consider this 
regime for the reasons given at the end of section 
II. 

The above pair of very narrow and very high 
boxes that touch each other produces an effec- 
tive potential which can be approximated very 
well by that resulting from a (^'-function inho- 
mogeneity. This can be shown by considering 
in Eq. (|T8|l the limit 6 — >■ 0, — > cx) with finite 
e := {2hY -h for a box pair at xi and X2 = xi + 2b 
in cell n = 0. The limit yields 



4e , X- 
— tanh J 



Xi 2^ 
— sech 



Xi 



I 



f/(°)(^,0 

(27) 

which is identical with U^^^ stemming from a 6'- 
function inhomogeneity gi = e6'{x — Xi). 



VI. COMPARISON WITH SIMULATIONS 

In order to check the predictions of the above 
collective coordinate approach with the opti- 
mized effective potential constructed in the pre- 
vious section, we have performed simulations for 
the full sine-Gordon system with the box 
inhomogeneities that produce the above effec- 
tive potential: two narrow and high boxes per 
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FIG. 4: Mean velocity vs. amplitude A of the ac force, (a) 
Numerical solution of the CC-Eqs. 151111 with U{X, I) given 
by Eg. 11811 . (b) Simulations of Eq.Q with the expression for 
V{x) given by Eqs.J^J and 1161 . V{x) is shown in the inset. 
The parameters used are l3 — 1, uj — 0.015, h = 25, b = 0.1, 
xi = 0.8, X2 = 1.0 



cell which have opposite signs and touch each 
other. We have numerically solved Eq.(^ us- 
ing the Heun scheme jsj] which yields the same 
results as the Strauss- Vazquez scheme js^, but 
considerably faster; another advantage of the 
Heun scheme is the possibility to include thermal 
noise jl^]. The spatial and temporal integration 
steps were Ax = 0.025 and At = 0.005. The spa- 



tial interval for the simulations had a length of 
190 units. The simulations were performed tak- 
ing as initial condition a kink soliton Eq.(j3]) at 
rest, which we allow to evolve. Then, after some 
transient time, we compute the mean velocity 
integrating over a period of time. This process 
is repeated again by varying the amplitude of 
the force, thus sweeping over the whole interval 
for the force amplitude. We have checked that 
the computation of the mean velocity does not 
change over several periods. 

The results in Fig^ show a very good qualita- 
tive agreement with the CC approach results in 
FigEt^- Interestingly the simulation results are 
more sensitive to the choice of the cell length 
than the CC approach results: both the max- 
imum Vmax of the average velocity and its po- 
sition change more when L is increased from 4 
to 8. The best value is Vmax = 0.172, which is 
350% higher than the best result 0.0382 for the 
6 -function inhomogeneities 0] . In this case the 
highest value of the soliton speed |f (t)| is about 
0.5 which has to be compared with the critical 
velocity c = 1. 

So far we have only considered the case of a low 
frequency {u = 0.015) of the driving force. Here 
both the (5-function and the box inhomogeneities 
produce a staircase in {v) vs. A. However, 
for higher frequencies the situation changes: For 
u = 0.05 the boxes still yield a staircase (Figl^I) 
while the (5-functions yield several windows ( see 
case (b) in section II). For u = 0.1 (not shown) a 
2-step staircase is seen for L = 4, while two and 
one windows are found for L = 6 and L = 8, 
respectively. 



VII. SUMMARY AND CONCLUSIONS 

We have studied how ratchet transport in inho- 
mogeneous sine-Gordon systems can be strongly 
enhanced, proceeding in three steps: 

1) In a collective coordinate approach the soli- 
ton dynamics was represented by a relativistic 
particle ratchet with an effective potential U, 
which is periodic and asymmetric. Using an 
ansatz in terms of Jacobi elliptic functions, the 
optimal U for the particle ratchet was found. 
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FIG. 5: Mean velocity vs. amplitude A of the ac force, (a) 
Numerical solution of the CC-equations O- Hllll . (b) Simula- 
tions of Eq.Q with the expression for V{x) given by Eqs.Q 
and llbH . The parameters used are /3 — 1, ui = 0.05, b — 0.1, 
h = 25, XI = 0.8, X2 = 1.0. 



a periodic array in which each cell contains two 
very narrow and high boxes with opposite signs 
that touch each other (this structure can be well 
approximated by a 6'- function). The optimiza- 
tion of the parameters (width and height of the 
boxes and the cell length) yielded high values for 
Vmax, the maximum of the average soliton veloc- 
ity (v) as a function of the driving amplitude 
A. For the frequency range [0.015, 0.1] Vmax is in 
the order of 0.15, where the soliton speed \v(t)\ 
reaches values which are not far away from the 
critical velocity c = 1. For frequencies higher 
than uj = 0.2 the ratchet effect vanishes. 

3) Simulations for the sine-Gordon system 
with the above inhomogeneities produce results 
for (v) vs. A which show a very good qualita- 
tive agreement with the results of the collective 
coordinate approach. 

From the results which we have obtained, 
we propose the following experiments on long 
Josephson junctions: Narrow inhomogeneities 
of opposite signs can be built by microresistors 
(critical Josephson current Jc decreased) and mi- 
croshorts (Jc increased). So far only equidistant 
arrays of either type of inhomogeneities were 
used, but we do not see problems in producing 
arrays in which microresistors are placed close to 
microshorts. 

A further enhancement of Vmax could be 
achieved by a replacement of the sinusoidal 
driving force f{t) by a more general periodic 
force (for homogeneous sine- Gordon systems an 
ansatz in terms of Jacobi elliptic functions, see 
below Eq. 25, yielded Vmax = 0.16 and 0.17 
in simulations and a CC-approach, respectively 
[s^l). Work on the joint optimization of the inho- 
mogeneities and the driving force is in progress. 



Uopt exhibits essentially two features: a) dou- 
ble peaks consisting of a positive and a negative 
peak with a steep slope in between, b) long flat 
parts between the double peaks. 

2) The inhomogeneities in the sine- Gordon sys- 
tem were chosen such that they produce an effec- 
tive potential which exhibits the above features 
of Uopt- It turned out that a very good choice is 
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